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LLR analysis at IfE - program LUNAR 


Goes back to mid of 1980ies (FESG, Munich) 

■ 2 main development phases: -1985-2001 (Munich), since 2006 at IfE 

■ LLR-only analysis 

- ephemeris computation based on initial values from DE421 (Sun, planets, 
largest asteroids) 

no further planetary data (optical, radar,...) included 

- ephemeris model: EIH equations of motion for all bodies as point masses 
Earth orientation: IERS conventions 

lunar orientation: Euler equations integrated together with translation 
including relativistic corrections (geodetic precession, Lense-Thirring) 

- additional forces due to inhomogeneous gravity fields, tides 

■ Combined analysis of solar system data and LLR planned in future 
(project in research unit FOR 1503) 









Program updates at IfE 


■ Data reduction compared with IMCCE, good agreement 

■ Results show some room for improvements -> updates in ephemeris 
computation needed 

- many modeling parts still from first development phase 

- simplifications (due to computation time reasons and accuracy 
requirements) 

- slight inconsistencies in force model (.interactions 1 ) resolved 

3 parts: Earth-related, additional gravitational effects, Moon-related 












Program updates at IfE - Earth related 


Earth tides - tidal acceleration 

- former model 

• single lag angle of time delayed tides 

• only Moon as tide generating body 

• effect on lunar translation 
new model (according to DE430 ephemeris) 

• degree 2 tidal potential (-98% of tidal effect) 

• arbitrary tide generating body possible 

• 5 tidal time delays (2 estimated, 3 fix at DE-values) 

• effect on Moon: 

■ via change in Earth's degree 2 potential coefficients on lunar 
translation and rotation 

■ via acceleration on lunar translation 
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Program updates at IfE - Earth related 


Earth tides - tidal acceleration 
- some results 

• Sun + Moon as tide generating body 

• tides from Jupiter, Venus -> <0.1 mm in r EM in 45 years 

• estimated time delays 

■ diurnal t 21 =575 s (DE430: 636 s) 

■ semi-diurnal t 22 =226 s (DE430: 219 s) 




Secular trend in C 20 

several models tested (linear, quadratic) 
best result in LLR analysis with model 
from IERS Conventions 2010 with 

C 20 =2.6xl0 _11 yr _1 



Year 


Roy and Peltier (2011) 
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Program updates at IfE - additional gravitational effects 


Interaction of Sun/planets with Moon 

- former model - interaction of 

• planets with point-mass Moon 

• Sun with lunar degree 2 

new model 

• planets with lunar degree 2 

• Sun with lunar degree 2+3 
main effect on lunar rotation, e.g. 

• ~ 19 mm on surface from Venus 

• ~ 4 mm on surface from Jupiter 
(from ephemeris with equal initial conditions) 

- in analysis: residuals decrease ~0.1 mm on average 
-> small effect but maybe needed in future 



EM 
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Program updates at IfE - additional gravitational effects 


Figure-figure interaction between Earth and Moon 

former model 

• simplified version of degree 2 - degree 2 coupling 

• effect on lunar rotation considered 

new model 

• coupling up to any degree/order of the 
gravitational field of Earth and Moon possible 
(Ilk, 1983) 

• effect on translation and rotation 




results 

• improvement due to complete 
degree 2-degree 2 coupling 

• Earth degree 2 - Moon degree 3 
-> some mm on surface 



Jahr 
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Program updates at IfE - Moon related 


Rotation of deformable Moon 

- former model 

• tidally deformed tensor of inertia only in rotation 

• no changes of potential coefficients (interaction with Earth, additional 
effect on translation and rotation) 

• no consistent core-implemention 
new model (according to 2-layer model of DE430) 

• basis tide-free tensor of ineria, elements from 

■ C 20 (GRAIL) 

■ C 22 , dynamical p (estimated in analysis) 

■ dynamical y derived 

• core moments with DE-fixed values for inertia-ratio C core /C Moon 
and core flattening 
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Program updates at IfE - Moon related 



• tidal + spin deformations on mantle tensor of inertia 

■ Love number k 2 fixed on GRAIL-value 

■ tidal deformation due to Earth (much larger than from Sun) 

■ 1 time delay used 

-> equations for rotation complete 
-> coupled differential equations for core + mantle 
-> coupling constant, initial rotation vector of core estimated 


r 


EM 



• degree-2 changes in potential coefficients from deformed mantle tensor 
-> enters in computation of external forces (translation/rotation) 
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Program updates at IfE - Moon related 


Effect of new rotation¬ 
modelling on post-fit 1 

residuals, wrms (without < 

core) 



Effect including core 
(which leads to improved 
libration modelling) on 

2-way residuals 
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Results - wrms 


1 -way wrms, comparison with former model 



Further investigations 

effects in longitude libration visible -^empirical correction as in DE430 
core flattening not yet estimated possible effect on librations 
- tidal and lunar rotation modelling room for improvement 
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Results - coefficients lunisolar nutation 


Periode MHB2000 [mas] old model new model (two test cases) 
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to larger discrepancies 
in other coefficients 
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Results - relativistic parameters 


Temporal variation of gravitational constant 


modelled as G(t) = G 
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moderate to strong correlations with 

• lunar core rotation vector (fixed in solution) 

• some station coordinates (constrained a little bit) 

- as single parameters: — = (0.7±0.8)xicr 13 yr -1 

G o 

— = (1.6 + 2.0)xl0 _15 yr _2 

G o 

estimated together: — = (0.8+l.l)xl0 _13 yr' 1 — = (-0.3±2.4)xio 15 yr -2 
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Results - relativistic parameters 


Equivalence principle 

estimating ratio A(m /m ) 

V 9 1 /EM 

• partials: computed numerically by introducing 
additional acceleration of Moon towards Sun 

• Afm/rn) = (-3.0±6.6)xl(T 14 

V 9 1 /EM 

Estimating Nordtvedt parameter rj 

• partials: analytical from synodic range 
variation: 13.1m*cos(D)* rj 

• ri = (-0.2 ± 1.2) x 10” 4 



into direction of galactic center (e.g. due to dark matter) 

• amplitude: A nc =3.0±3.3mm 

• additional acceleration a =(-i.i+i.2)xio -6 
in parts of 1.9e-8 cm/s 2 9C 
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Results - relativistic parameters 


■ PPN parameters p, y 

- included in EIH-equations of motion, partials numerically 

- correlated with station coordinates (constrained) 

show also correlations with z-coordinate of lunar initial orbit values 

- p-l = (0.9±1.0)xio -4 

- y _i = (_i 2 + 1 . 6 )xicr* (not as accurate as Cassini-result) 

p from combination of PPN-parameters 0.25(y+r|+3) and Cassini- y 
P-l = (0.03±6.1)xl0“ 5 

■ PPN preferred frame a,, a 2 w.r.t. cosmic microwave background 

- a, =(-l.l±2.0)xl(T 5 

a = (-0.6±0.9)xl0 -5 ( n °t as accurate as test with Sun's spin) 









Results - relativistic parameters 


Geodetic precession of lunar orbit 

- introducing GP a second time as additional acceleration 

factor h gives relative deviation in from Einstein's theory (~1.9 as/cy) 
strong correlation with 

• lunar core rotation vector (fixed) 

• z-component of lunar initial velocity (fixed) 

h = (-0.6 + 2.0) x 1CT 3 

■ Yukawa-term (1/r 2 -test), acceleration due to v^ M — — G MeMm (i+ae~ r/1 ) 

- interacting range X = 380000km 

- coupling constant a estimated 
correlations and fixed values like GP 

- a = (-4.0 ± 5.0) x 10” 12 










Summary 


IfE-LLR ephemeris model updated 

- tidal acceleration, secular trend in Earth's C 20 

- additional gravitational interactions planets-Moon, Earth-Moon 
lunar rotation as 2-layer core/mantle model 

- effect of lunar deformation in all lunar equations of motion (translation 
plus rotation) 

■ Smaller residuals and more accurate parameter estimation 

- increased accuracy in relativistic parameters (strong limits for validity of 
equivalence principle and gravitational constant) 

no deviation from Einstein's theory of gravity up to now 








